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ABSTRACT 



The use of differentially deflected fairwater planes 
to control submarine rolling is studied. Because of coupling 
between pitch and roll angles, the snap roll that occurs in 
a high speed hard turn affects the stability of a submarine 
not only in the horizontal plane but also in the vertical 
plane. Direct roll control was achieved by making use of the 
fairwater planes in a differentially deflected mode such that 
they could give counter moment to reduce the snap roll. A 
roll controller was designed as a position and velocity feed- 
back controller. Controlled roll angle improved the depth 
and pitch stability of submarine. 
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I. INTRODUCTION 



The importance of automatic control of oitch-depth and 
roll of a turning submarine has become more obvious with the 
improvements of high speed nuclear submarines. This study was 
aimed at smoothing the roll of a turning submarine in hioh 
speed, which is known as snap roll and causes the basic problem 
Because of coupling between pitch and roll, before attemoting 
to smooth the snap roll, controlling of pitch and depth in 
high speeds turns was also studied. In this study an auto- 
matic course controller was not considered. 

The methods and techniques, which can be used to reduce 
and smooth the snap roll, can be categorized under two head- 
ings. The first concerns changes in the naval architectural 
characteristics of the designs, such as increasing GM (meta- 
centric height) and reducing sail size. The second category 
involves those alternatives which make use of an automatic ship 
control system. Rudder sequencing and speed reduction fall 
under this classification. Reference 1 investigates the affect 
of increasing GM, reducing sail size and speed. In Reference 2 
Stamps designed on automatic roll controller which makes use of 
rudder sequencing as a function of aooroach speed and instan- 
taneous roll angles. 

In recent studies, which were made by Naval Ship Research 
and Development Center (NSRDC) , it was proposed to enhance the 
control of roll in high soeed turns by using the Faiir^vater 
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planes differentially deflected so that they can be used to 
give a counter moment to reduce the roll. The investigation 
of the differentially deflected Fairwater planes effect was 
part of the project called "Improved Control For Advanced 
Submarines". The project was carried out under Program Element 
No. 62754N and Task Area 2F434001. The work unit number was 
1-1563-001-74. In Reference 3, estimating the effectiveness 
of differentially deflected sailplanes was investigated as a 
part of project mentioned above. Reference 3 was the inspira- 
tion of this thesis and using differentially deflected sail- 
planes for direct control of roll was chosen as the design goal. 
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I 



II. THE NATURE OF THE PROBLEM 



Turning characteristics of a surfaced submarine, more 
or less, looks like those of a surface ship. But the situation 
in the submerged position shows big differences. These dif- 
ferences are the result of different naval architectural char- 
actacteristics . Sail structure can be considered the main dif- 
ferences and the main source of roll problems. It is v/ell 
knov/n that when a surface ship goes into a turn, it experiences 
an initial inboard roll. After a very short transient it heels 
outboard. The reversing of the roll anqle is primarily due to 
an asymmetric rudder. Since the submarine rudder is generally 
well submerged, the surfaced submarine has an outboard roll 
angle during both the initial and the steady state phase of a 
turn. For a typical submarine, at moderate speed this outboard 
roll angle is less than 10° and does not cause any big problem. 

If a submerged submarine goes into a high speed turn, 
differences from the surfaced behavior are noticed and the 
problem becomes three-dimensionalized. In the very first phase 
of the turn, the ship has a small lateral velocity v and a 
small rotational velocity r. For a starboard turn, at some 
point along the centerline, XI, the lateral velocity due to 
r, rXl, equals the lateral velocity - v. (For symbols see 
Figure 1, 2, and Appendix A). The point, at which these two 
velocities are equal to each other, is called the instant 
center of turn. At some point, forward of the instant center, 
X 2 , the lateral velocity, -v, gives a velocity component from 
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the starboard side. Along with the ahead velocity coinoonent 
of the ship, these velocities contribute an angle of attack 
from starboard. At this point, if the submarine has an append- 
age (sail) , a lift force is produced on the sail directed to 
port. Depending on the sail configuration and the large moment 
arm of the sail. This force, directed to the port gives a small 
outboard (port) roll. In many submarines this outboard roll can 
not be felt. 

In the steady state phase the angle of attack shifts to 
the other side (port) and results in an inboard roll. But con- 
trary to the outboard roll this inboard roll is quite signifi- 
cant and at high speeds it can exceed 30° which is considered 
hazardous to both men and equipments. This inward roll is 
called the SNAP roll. In Figure 5, the roll characteristics 
of the ship to a constant 35^ left rudder angle at 24 knots is 
shown. As can be seen, the initial outboard roll (in the 
simulation, outboard is the starboard side) is very small, 
less than 1°. But the snap roll reaches 37°. 

Snap roll is not the only problem that can be faced when 
the submarine goes into a turn. Because of the appendage 
(which is mainly the sail) , in the submerged position, sub- 
marines do not have very good hydrodynamics from this body 
structure. When the ship experiences a submerged turn the 
sail mainly destroys the waterflow around the hull. As a 
result a pressure difference is created between the upper 
Dortion of the hull and the bottom portion of the hull. 
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And a normal force results directed dov/nward which bodily 
pushes the ship down. It is recalled that, in a turn, after 
the initial phase, the submarine has an inward roll with the 
deflected rudder. Rolling inward, causes this deflected 
rudder to act like a sternplane which gives the submarine a 
down pitch movement. All of these effects come together and 
cause the ship to rise or dive depending on the peculiarity 
of its design. As it is stated in Reference 4, "The exact 
mechanism of this is not well understood". In Figures 3, 4, 
and 6, the depth, oitch, and rudder response of the ship, to 
a 35° left rudder angle at 24 knots is shown. It is interest- 
ing to indicate here that after a long transient period, the 
submarine starts to rise with the bow down eventhough it dives 
initially. It should be remembered that in this simulation of 
the model, no control surface (stern and fairwater plane) was 
used. The changes of these characteristics will be seen after 
the depth-pitch and roll controller is designed. 

The complexity of controlling the turning submarine 
mainly comes from the coupling between the roll and pitch. 
Since roll and oitch are coupled, to be able to control the 
roll, which is our primary concern, pitch control to some ex- 
tent must be accomplished. In the following sections a pitch 
and depth controller is designed by using only the sternolanes 
because of the goal of using the fairv/ater plane as a main 
part of the direct roll controller. 
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Figure 1. Hydrodynamic Forces In A Turn. 
(After Reference 4) 
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Figure 2. Submarine Axes 
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Figure 3. Depth vs. Time. With No Controller 

UCK = 24 Knots. Rudder Ordered = 35°. 
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Figure 4. 



Pitch vs. Time. With No Controller. 
UCK = 24 Knots. Rudder Ordered = 35°. 



23 




Figure 5. Roll vs. Time. With No Controller. 

UCK = 24 Knots. Rudder Ordered = 35°. 
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Figure 6. Rudder Response vs. Time. With Mo Controller. 
UCK = 24 Knots. Rudder Ordered = 35°. 
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III. DIGITAL COMPUTER SIMULATION OF SUBMARINE MOTION 



The mathematical model of the totally submerged sub- 
marine is represented by the equations of motion x^^hich consist 
of six equations such that each of them represent the force 
and moment equations on and around the three rectangular co- 
ordinate axes. These equations were derived by NSRDC in 
Reference 5 and for convenience was repeated in Apoendix A. 

In this thesis, derivation of these equations was not our 
concern. More pronounced knowledge about that can be found in 
References 4, 5, and 6. 

The equations of motion, which represent the motion of 
the totally submerged submarine, in six degree of freedom, con- 
sist of Hydrodynamic coefficient, in undimensionaliced form 
except the m (mass) and I , I , I (three moments of inertia). 

X V Z 

In Reference 7, Drurey translated these six equations into the 
DSL (digital simulator language) form and originated the com- 
puter program which is the mathematical model of the ship. 

For the necessity of using the hydrodymanic coefficients in 
undimensionalized form, the equations containing I , I , I , 

X ^ 2 

both sides are divided by 1^ and these three moments of 
inertia were used in undimensionalized form. For the same 
reason the equations containing m were divided by 1^. 

Since the trim control of the ship was not considered 
in this study the ship was assumed in the trimmed condition. 
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For trimming the ship the hydrodynamic coefficient 
Z* and M* were set to zero. These two hydrodynamic co- 
efficients represent the force and moment acting on the ship 
when the control surfaces are at zero deflection. 

Since all of the simulations of this study were done 
in a turn, an accurate simulation must include the rudder 
actuator dynamics such that any phase lag could not produce 
any unstability. Actuator dynamics used in this thesis v;as 
the same as used in Reference 3 with the rudder deflection 
angle 5°/sec., so that at the end of the study it could give 
acceptable comparison level between the two designs. A block 
diagram of the rudder actuator is shown in Figure 7. 

Since our prime concern is to control the roll as v/ell 
as the depth and the pitch when the submarine is in a turn, a 
series of tests were run at different speeds to validate the 
models originated in Reference 7 by Drurey. The model was 
highly satisfactory as shown in Figures 8 through 25. These 
figures record responses in depth, pitch, roll, yaw, speed 
changes, and rudder angle at 18, 12, and 6 knots. In this 
test, turn was commanded 10 seconds after the simulation 
started, to validate the initial condition response of the 
model. Since the submarine was trimmed by eliminating Z* and 
, it is seen that until the turn was commanded all responses 
are zero. Simulation at 24 knots was already shown in Figures 
3 through 6 . 

It was observed that, at 6 knots the ship continues 
diving unlike the simulation of 12, 18, and 24 knots. The 
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reason for this was that at low speed enough pressure dif- 
ference to cause the ship to rise slov/ly could not be created. 

After the validation of the model it was decided to pro- 
ceed with the designing of the controllers. In the following 
sections first, the depth and pitch controller design was stud- 
ied. An optimal control scheme such as in Reference 7 was used. 
The sternplanes was the control surface of the controller, and 
the fairwater plane was saved for use in the roll controller. 
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Figure 7. Rudder Actuator Model. (After Reference 2) 
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Figure 8. 



Depth vs. Time. With No Controller. 
UCK = 18 Knots. Rudder Ordered = 35°. 
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Pitch vs. Time. With No Controller. 
UCK = 18 Knots. Rudder Ordered = 35*^. 
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Figure 10. Roll vs. Time. With No Controller. 

UCK = 18 Knots. Rudder Ordered 35*^. 
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Figure 11. Speed Change vs. Time. With No Controller. 

UCK = 18 Knots. Rudder Ordered = 35°. 
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Figure 12. Yaw vs. Time. With No Controller. 

UCK = 18 Knots. Rudder Ordered = 35°. 
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Figure 13. Rudder Response vs. Time. 

UCK = 18 Knots. Rudder Ordered = 35*^. 
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Figure 14, 



Depth vs. Time. With No Controller. 
UCK = 12 Knots. Rudder Ordered = 35°. 
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Figure 15. 



Pitch vs. Time. With No Controller. 
UCK = 12 Knots. Rudder Ordered = 35°. 
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Figure 16. Roll vs. Time. With Mo Controller. 

UCK = 12 Knots. Rudder Ordered = 35°. 
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Figure 17. 



Speed Change vs. Time. With No Controller. 
UCK = 12 Knots. Pudder Ordered = 35°. 
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Figure 18. Yaw vs. Time. With No Controller. 

UCK = 12 Knots. Rudder Ordered = 35*^. 
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Figure 19. Rudder Response vs. Time. With No Controller. 
UCK = 12 Knots. Rudder Ordered =» 35^. 
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Figure 20. Depth vs. Time. With No Controller. 

UCK = 6 Knots. Rudder Ordered = 35°. 
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Figure 21. 



Pitch vs. Time. V7ith No Controller. 
UCK = 6 Knots. Rudder Ordered = 35°. 
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Figure 22. 



P.oll vs. Time. With No Controller. 
UCK = 6 Knots. Rudder Ordered = 35°. 
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Figure 23. Speed Change vs. Time. With No Controller. 
UCK = 6 Knots. Rudder Ordered = 35°. 
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Figure 24. 



Yav; vs. Time. With No Controller. 

UCK = 6 Knots. Rudder Ordered = 35°. 
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Figure 25. Rudder Response vs. Time. With No Controller. 
UCK = 6 Knots. Rudder Ordered = 35°. 
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IV. AUTOflATIC DEPTH AND PITCH CONTROL 



A. OPTIMAL SOLUTION TO THE LINEARIZED SUBMARINE EQUATIONS 

The necessity of controlling the depth and pitch of a 
turning submarine, as was stated before, stems from the coupling 
between the states such that any changes of one of them directly 
affects the other one. For this reason, depth and pitch control 
to some extent must be accomplished so that it can provide a 
good stage for the roll control. There are many ways of con- 
trolling the depth and pitch which depend upon the way of using 
stern and fairwater control surfaces combination. The scheme 
that was used in this thesis was an optimal control way which 
was originated by Drureys in Reference 7. To preclude unneces- 
sary repetition, only guidelines of the method and the differ- 
ences of this study are to be discussed here. More pronounced 
knowledge can be found at Reference 7 and 8. 

It can be seen in Appendix A that the equations of motion 
are nonlinear. To be able to use this set of equations, they 
must be linearized. This linearization was done in Reference 7, 
and with one control input (sternplane) they can be written in 
the form of: _ 
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If they are put in state variable form 
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— 




— — 


w 


=r 


^11 ^21 




W 




®11 


• 

Q 




^12 ^22 


1 


Q 




®12 



DS 



Where 



Determ = jJ(Iy-M<5) (m-Z^) - 
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[>1«,zq + <'”-Z«'>Mq] “/I. Deterra 

j^(I^-M(5) Zds + Z^MdsJ ^ /\ * 
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The problem was thought as a linear regulator problem. The 
general scheme of linear regulator problem is shown in Figure 
26. 

The cost function is 

j = * I Te^'^QE + u'^Ru"] dt 



= j ( + u'^RuJ 



E = X-r and r is the reference vector which is 



ordered depth 
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ordered pitch 
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Figure 26. Plant And Optimal Feedback Controller For Linear Regulator Problem. 
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As a result 



E = X can be written. 
Following the Reference 7 and 8 



K = -KA+KBR“^b'^K-Q-a’^K and 
— 1 T 

U = -R B KE can be written. 

If the linear equations are augmented and put in the state 
equation form with the definition 



E = X = 
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In this study, in order to make the depth rate a state vari- 
able the following transformation was used 

Depth rate = w-u* 0 where v! represent the comnonent of 
u in the z direction, Q represents the pitch angle. In 
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Reference 7 this transformation was originally 

Depth Rate = w and the differences between these two 

transformations showed itself in the calculation of K values 
which is discussed in the following cages. 

Referring to the cost function 
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The selection of the weighting factors was a trial and 

error process. The series of weighting factors depending upon 

the relative severity of the influence of the states to each 

other was tested and as a result 

E = Depth Error Weighting = 10 

D = Pitch Error Weighting = 8000 

C = Control Incut Weighting = 100 

A = B = 0 was chosen. 

After solution of n(n+l) > differential eauations , 

'2 

K values associated vrith the feedback aain was found. It is 
observed that steadstate values of X are constant which is a 
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convenient condition for optimality. Results of this solution 
are shown in Figures 27 through 36. For comparison the original 
A matrix of Reference 7 with the transformation Depth rate = w 
was also used to calculate K values. One of the results is 
shown in Figure 37. It is observed that it could not reach 
steady state value in the reasonable time period. 

After K values associated with the feedback gains were 
found, necessary gains were found via the optimal low. 

1 T 

u = -i B K E 

The result was 

u = DSAD = -0. 31623*Z0ER-1. 792*Z0 d0T + 

36.069*PERR = 102. 63*P1D0T. 

All of the feedback gains were found at 15 knots speed. 

It is obvious that these gains associated with the hydro- 
dynamic coefficient used are function of the speed. To make 
the gains compatible with the soeed range of the ship, they 
must be scaled with the function of the speed. In Reference 2, 
the controller gains adjustment as a function of the speed was 
discussed deeply and found that gains associated with the depth 
rate and pitch role error channel was inverselv proportional to 
the instantaneous speed as gains associated with the depth 
and pitch error channel remains same. After the justification 
of these results 

DSAD = -0.31623*ZOER-1.792*(15) / *ZoDOT+ 

'^UK 

36.069*PERR + 102.63*(15)/ *P1D0T 
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27. Kll vs. Time. UCK = 15 Knots. 

D = 8000, A = 0, E = 10, B = 0, C = 100. 
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Figure 28, K12 vs. Time. UCK = 15 Knots. 

D = 8000, A = 0, E = 10, B = 0, C = 100, 
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Figure 29, Kl3 vs. Time. UCK = 15 Knots. 

D = 8000, A = 0, E = 10, B = 0, C = 100. 
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Figure 30. Kl4 vs. Time. UCK = 15 Knots. 

D = 8000, A = 0, E = 10, B = 0, C = 100. 
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Figure 31. K22 vs. Time. UCK = 15 Knots. 

D = 8000, A = 0, E = 10, B = 0, C = 100. 
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Figure 32. 



K23 vs. Time. UCK 
D = 8000, A = 0, E 



15 Knots. 

10, B = 0, C = 100. 
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Figure 33. K24 vs. Time. UCK = 15 Knots. 

D = 8000, A = 0, E = 10, B = 0, C = 100. 
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Figure 34. K33 vs. Time. UCK = 15 Knots. 

D = 8000, A = 0, E = 10, B = 0, C = 100. 
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35. K34 vs. Time. UCK = 15 Knots. 

D = 8000, A = 0, E = 10, B = 0, C = 100. 
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Figure 36. K44 vs. Time. UCK = 15 Knots. 

D = 8000, A = 0, E = 10, B = 0, C = 100. 
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Figure 37. Kll vs. Time. UCK = 15 Knots. 

With Transformation Of Depth Rate = W. * 
D = 8000, A = 0, E = 10, B = 0, C = 100. 
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*Z0D0T + 



= -0. 31623*Z0ER-26. 38 . 

^UK 

36.069*PERR + 1539.45 , *P1D0T 
where UK is the instantaneous soeed (knots). 

B. AUTOMATIC DEPTH-PITCH CONTROLLER AND SIMULATION OF THE 

RESULTS 

The feedback gains, which are the ootimal solution of 
the linearized vertical plane equations must be put into the 
general controller schema in which the submarine dynamics are 
represented in six degrees of freedom (i.e., non-linear 
equations of motion) . Because only the simulation of this can 
justify the validity of the solution. The complete modified 
version of the Drurey's depth and pitch controller is shown in 
Fighre 38. As it is shown the sternplane controller (actuator) 
was not designed as a part of the controller and put into the 
controller separately with the plane rate 7°/sec. The dynamics 
of the actuator is the same as Drurey's and Stamps' controller 
used. The same actuator dynamics also was used in the roll 
controller which is to be discussed in the following section. 
The actuator block diagram is shown in Figure 39. To preclude 
stability which occurs as a result of the excessive uses of 
the sternplane and avoid use of the control surface with big 
deflections the limiter was put into the pitch and depth error 
channels. Pitch error limiter was 10° and depth error limiter 
was 20 feet and found by a trial and error process. 

After the completion of the controller design, various 
tests were run to see the effect of the pitch and depth con- 
troller on depth, pitch, roll and sternplane responses. 
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Figure 38. Depth And Pitch Controller. 
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(DS Refers To Stern Planes Actuator. DF Refers to 
Fairwater Planes Actuator) . 



B 





The results of tests at 6, 12, 18, and 24 knots to a constant 
35° left rudder angle are shown in Figures 40 through 55. 

^^hen these results are comoared with the responses of the ship 
with no pitch and depth controller which is shown in Figures 
3 through 6 at 24 knots and in Figures 8 through 25 at 18, 12, 
and 6 knots. The following comoarisons can be drawn; 

1. At 12 and 18 knots controller kept the ship's depth 
and pitch stable after very reasonable transient tiipe. 

The steady state pitch and depth are (5.2)° - (9.2) feet 
at 12 knots and (4.7)° - (9) feet at 18 knots. 

2. At 6 knots, the controller failed to keep the pitch and 
depth stable. The shio very slowly kept on losing depth. 
The reason was insufficient speed scaling of the feedback 
gains and needs for another control surface (Fairwater 
planes) . Lack of fairwater plane control surface in the 
depth-pitch controller unabled it to keep ship stable at 
low speeds. But in this thesis the primary concern was 
to minimize the snap roll which occurs at high speed. 

In Reference 2 Stamps collected data which shows the peak 
snap roll as a function of the approach speed and ordered 
rudder angle. This data is repeated at the following 
page for convenience (Table I). It is seen that dangerous 
snap roll starts occuring at approach speeds above 12 
knots. For this reason unability of the controller at 
low speed was neglected. Because, a control scheme, which 
at low speed uses Drury's original depth-pitch controller 
(Fairwater plane is part of controller) and at high 
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speed uses this design (’^airwater plane is not a Dart 
of depth-pitch controller but is a part of the roll con- 
troller) can be designed and switching from one to another 
can be achieved. 

3. At 24 knots the controller again failed to keep the ship 
stable. Depth and pitch response went into oscillation 
with big amplitudes. Original roll response of the ship 
which is shown in Figure 5 was destroyed in the sense of 
decreasing amplitudes of roll oscillation. It did not 
reach stable value and appeared to be oscillating in the 
range of 10°. The reason was the following. The big snap 
roll which was around 37°, initially gave very big dis- 
turbance and oscillation. With these big disturbances 
the controller which uses only the sternolane as a con- 
trol surface was unable to do the job. It was thought 
that, if there had been any controller which could have 
decreased the snap roll (Roll controller) it would have 
been able to stabilize the ship in pitch and deoth as 
well as in the roll response. 

Based on the results of these tests, proceeding with the 
roll controller design which would make use of fairwater plane 
as a control surface was decided. 
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TABLE I 



PEAJ< SNAP ROLL ANGLE (DEGREES) VS. APPROACH SPEED (KNOTS) 





AS A FUNCTION 


OF ORDERED 


RUDDER 


ANGLE (DEGREE) 




SPEED 

RUDDER 


(KNOTS) 4 


8 


12 


16 


20 


24 


1 


0.15 


0.66 


1.45 


2.58 


4.03 


5.83 


2 


0.24 


0.99 


2.22 


3. 95 


6.18 


8.88 


3 


0.31 


1.24 


2.80 


5.00 


7. 84 


11. 28 


4 


0.37 


1.45 


3.29 


5.87 


9.26 


13.27 


5 


0.41 


1. 63 


3.71 


6.71 


10.32 


15. 72 


6 


0.45 


1. 80 


4.04 


7. 39 


11.65 


17.90 


7 


0.49 


1.93 


4.46 


7.93 


13.11 


19.87 


8 


0.52 


2. 09 


4.81 ■ 


8.36 


14.47 


21.64 


9 


• 0.55 


2.22 


5.12 


9.10 


15.72 


23.23 


10 


0.58 


2.33 


5.38 


9.88 


16.88 


24. 65 


11 


0.60 


2.42 


5.60 


10.62 


17.96 


25.94 


12 


0.63 


2.50 


5. 78 


11. 31 


18.95 


27.13 


13 


0.65 


2.60 


5.94 


11.97 


19.86 


28.20 


14 


0.67 


2.72 


6.06 


12.58 


20.71 


24.17 


15 


0.68 


2.83 


6.33 


13. 17 


21.49 


30.04 


16 


0.70 


2.92 


6.65 


13.72 


22.20 


30.83 


17 


0.72 


3.01 


6.96 


14.23 


22.87 


31.57 


18 


0. 74 


3.09 


7.25 


14.72 


23.48 


32.22 


19 


0.76 


3.16 


7.54 


15.16 


24.04 


32.81 


20 


0.77 


3.22 


7. 80 


15. 59 


24.55 


33. 34 


21 


0.78 


3.27 


8.06 


15.98 


25.02 


33. 82 
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SPEED (KNOTS) 
RUDDER 


4 


8 


12 


16 


20 


24 


22 


0.79 


3.32 


8.30 


16.34 


25.45 


34.25 


23 


0.80 


3.36 


8.52 


16.68 


25.84 


34.64 


24 


0.81 


3.39 


8.74 


16.99 


26.20 


34.99 


25 


0.82 


3. 42 


8.95 


17.28 


26.51 


35.30 


26 


0.83 


3.44 


9. 14 


17.55 


26.80 


35.57 


27 


0.84 


3. 46 


9.32 


17.79 


27.06 


35.81 


28 


0.86 


3.47 


9.49 


18.01 


27.29 


36.02 


29 


0.87 


3.48 


9.65 


18.22 


27. 50 


36.20 


30 


0.87 


3.49 


9.80 


18.40 


27. 68 


36.36 


31 


0. 88 


3.49 


9.94 


18.57 


27. 84 


36.50 


32 


0.89 


3.51 


10.06 


18.71 


27.97 


36.60 


33 


0.89 


3.58 


10.18 


18.85 


28.09 


36.68 


34 


0.90 


3.64 


10.29 


• 18.96 


28.19 


36.75 


35 


0.90 


3. 70 


10.39 


19.06 


28.27 


36.80 
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Figure 40. Depth vs. Time. With Denth-Pitch Controller. 
UCK = 6 Knots. Rudder Ordered = 35°. 
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Figure 41. Pitch vs. Time. With Depth-Pitch Controller. 
UCK = 6 Knots. Rudder Ordered = 35*^. 
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Figure 42. 



Roll vs. Time. With Deoth-Pitch Controller. 
UCK = 6 Knots. Rudder Ordered = 35°. 
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Figure 43. Sternplane Angle vs. Time. ""Tith Depth-Pitch 
Controller. UCK = 6 Knots . Rudder Ordered 
= 35°. 
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Figure 44. Deoth vs. Tine. 

UCK = 12 Knots. 



With Death-Pitch Controller. 
Rudder Ordered = 35°. 
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Figure 45. 



Pitch vs. Time. With Depth-Pitch Controller. 
UCK = 16 Knots. Rudder Ordered = 35^. 
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Figure 46. Roll vs. Time. With Deoth-Pitch Controller. 
UCK = 12 Knots. Rudder Ordered = 35°. 
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Figure 47. Sternplane Ancle vs. Time. With Deoth-Pitch 
Controller. UCK = 12 Knots. Rudder . Ordered 
= 35°. 
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Figure 48. Depth vs. Time. With Deoth-Pitch Controller. 
UCK = 18 Knots. Rudder Ordered = 35°. 



80 



.00 4.00 6.00 8.00 10.00 






X SCALE = 40 Seconds Per Inch 
Y SCALE = 2 Degrees Per Inch 




Figure 49. 



Pitch vs. Time. With Depth-Pitch Controller. 
UCK = 18 Knots. Pudder Ordered = 35°. 
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With Depth-Pitch Controller, 
. Rudder Ordered = 35°. 



Figure 50 



Roll vs. Time. 
UCK = 18 Knots 
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Figure 51. Sternplane Angle vs. Time. With Depth-Pitch 
Controller. UCK = 18 Knots. Rudder Ordered 



= 35°. 
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Figure 52. Depth vs. Time. 

UCK =24 Knots. 



V7ith Depth-Pitch Controller. 
Rudder Ordered = 35^. 
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Figure 53. Pitch vs. Time. 

UCK =24 Knots. 



With Depth-Pitch Controller. 
Rudder Ordered = 35°. 
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Figure 54. Roll vs. Time. With Depth-Pitch Controller. 
UCK = 24 Knots. Rudder Ordered = 35°. 
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Figure 55. 
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vs. Time. With Deoth-Pitch 
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V. AUTOMATIC ROLL CONTROLLER 



In the introduction of the thesis, the possible choices 
of controlling roll was • discussed briefly and stated that 
fairwater planes v;ere to be used as control fins to control 
the roll. The reason that brought uo this idea lies behind 
the needs of more control surface than was available in the 
present design criteria. Roll control of a high speed sub- 
marine has continued to be a problem over the years because 
of limited control surface. Control engineers have been re- 
stricted using rudder, stern and fairwater plane. All of the 
control surface was not meant to control roll. In the present 
design criterias stern and fairwater planes are meant to con- 
trol the submarine in the vertical olane motion (depth-pitch) . 
The rudder has been used as a part of course controlling of 
the ships. Among the three control surfaces, rudder is the 
only one that has direct correlation with the submarines roll 
angles. In the latest study, which was done by Stamps 
/Reference ^7/ it was used to control roll. The concept of 
Stamps roll controller was based on the idea of adjusting the 
initial rudder angle order as a function of the integral of 
error between the allowed maximum roll and actual roll angles. 
The initial rudder order was chosen such that the peak roll 
expected for a given approach speed would be less than the 
maximum allowed roll. The integral of roll error was then 
computed and scaled to represent and additive term applied to 
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the initial rudder order. Stamps design didn't give any 
structural changes to the present navy submarines design 
criteria. In a sense of simplicity it was perfect. But it 
prohibited using hard and excessive rudder order which is 
highly desirable at some submarines required maneuvering in 
certain tactical areas. Since the Stamos design slowed down 
the yaw rate and caused the ship to change its course very 
slowly, alternative design ways were investigated and using 
fairwater plane in the differentially deflected mode was 
chosen as a possible improvement. 

The concept of roll control by means of fairwater planes 
is based on the idea of deflecting the fairwater plane dif- 
ferentially such that it can give a roll moment in opposite direc- 
tion to the instantaneous roll angle. If the ship has a 
roll angle to starboard side, the olanes are to be deflected 
to give a roll moment to portside. The positive sense of this 
additional roll moment created by the differentially deflected 
sailplanes is the same as of Reference 5. The positive sense 
of sailplane deflection angle adopted in this thesis (positive 
for port sailplane deflected leading edge up, starboard sail- 
plane deflected leading edge-down) is in agreement with that 
used for other control surfaces in Reference 5: Positive 

when the surfaces are deflected in such a direction as to in- 
crease the relevant angle of the submarine about its mass 
center. In the case of differentially deflected sailplanes, 
the angle is the roll angle and this is defined in Reference 5 
to be positive starboard side down. 
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Since, in the present design criteria of navy submarines, 
sailplanes are not used in a differentially deflected mode, 
standard equations of motion which are developed by NSRDC in 
Reference 5, do not have correlative terms which take the 
additional moment term, created by the sailplanes into account. 
For this reason, before starting to design the roll controller, 
this additive moment must be estimated. Once the counter 
moment is found it can be placed in the righthand side of the 
roll moment equation. 

A. ESTIMATING ROLL MOMENT DUE TO DIFFERENTIALLY DEFLECTED 

SAILPLANES 

In vertical plane motion, the normal force due to the 
fairwater planes (sailplanes) is given by the equation of 

Fjj = L^u^-Z-^ • 6b 

Where ^ = 2 

L = ship length in feet 
U = forward speed 

2 6b : hydrodynamic coefficient associated with the 
sailplane deflection (in conventional mode - right and left 
side moves together in the same direction) . 

6b = Deflection of the sailplane in radians. 

After setting =2 the normal force equation becomes 

Fj, = L^-U^-ZgySb 

This force is due to both sides of the sailplanes. The 
force due to one side of the plane pair is half of the total 
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force that is j Sb 

If both sides of the planes 

are deflected by the same amount and in opposite direction 
compared to each other, this configuration creates a moment 
associated with the moment arm between these opposite forces. 
This is shown in Figure 56. Total amount of the moment due to 
sailplanes can be written 

M 

M = L^* • Z^Iq* Sb. (moment arm) 

= L^‘U^*Sb» arm) 

if this momemt term is placed in the right side of the 
equation of motion about the body axis system X-Axis (roll 
axis) , the complete equation takes the form of 

lx? + (I^-lY)qr = /~Kp' p + Kqr'qr + Kr'r + P|P1_7 

+ /“Kp'up + Kr’ur + K^'v 

+ /“K*’u^ + K^’uv + K^j^j'v|(v^ + w^)| 

+ + L^-u2.K5^.Sr 

3 2 

+ L *U • (Z^^‘ moment arm 
+ Sin Cos 

In this last form 6b is replaced by to indicate that the 
fairwater plane is to be used in differentially deflected mode 
( ^^represents deflection anqle of fairwater plane in con- 
ventional usage. Both right and left sides move in the same 
direction). 



>A‘^F 



l-L“-U^-ZCj^ -Sb- (moment arm) 



= 2 * 
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Figure 56.1. 

. 2 . 



Positive Direction Of Differentially 
Deflected Sailplanes. 

Negative Direction Of Rolling Moment. 
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The moment arm used in this study was estimated by in- 
voking the data from Reference 6 that belongs to the un- 
classified fictitious submarine. In the reasonable proportion 
to the length of the ship, the moment arm was chosen as 7 ft. 
and, according the data, it was thought acceptable. 

If the new hydrodynamic coefficient is defined associated 
with the fairv^ater planes in differentially deflected node 




moment arm) 
L 



% = (0.00558 X 

Z^P = 0.0015517. 



After the estimation of Z^p, as a hydrodynamic co- 
efficient, it is placed in the modified equation discussed 
above. After these modifications, in the equations of motion 
all of the term which is a function of the Sb (sailplane de- 
flection in the conventional mode) was set to zero by defining 
Sb= 0. As result, the equation of motion about the X-axis 
(roll axis) has correlative term between the fair^^ater plane 
and roll angle of the ship so that this additive moment term 
can be used to give feedback to smooth the snap roll. 



B. ROLL CONTROLLER DESIGN 

In contrast to the depth and pitch controller, the roll 
controller was designed by using the nonlinear equations of 
motion and linearizing was not attempted. In the previous 
section, the modification of the roll moment equation due to 
the imposition of moment term which stems from the fairwater 
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planes deflection was discussed. All through the roll con- 
troller design this modified version of the equations was 
used. The reason that led to using six nonlinear equations 
was such that linearizing of this equations in six dereee of 
freedom was found difficult because of the terms coupled with 
each other and mainly was thought that from the outout of the 
system dynamics the roll angle and the roll rate (p) could 
give enough feedback information to accomplish the compensation. 

The principles for designing the roll controller were 
simple. If the ship has a positive roll angle (starboard) 
the fairwater plane is to be deflected so that it can give 
negative roll moment, i.e., starboard sailplane is to be de- 
flected leading edge up while port sailplane is deflected 
leading edge down or vice versa. As long as the ship has a 
roll angle this would cause the fairwater plane deflection by 
way of the feedback channel used. In the design it was assumed 
that, the faimi^ater planes actuater was capable of givina e- 
qual amounts of deflection command in opposite directions to 
both starboard and port sailplanes. Since the design was to 
invoke extensive simulation study (since nonlinear equations 
were used) in determination of the feedback parameters, the 
first basic controller attempt was the proportional controller 
because of its simplicity to design and implement. The pro- 
portional controller is shown in Figure 57.1. Referring to 
the figure, if the submarine has a roll angle this ^-/ould give 
the position feedback to force the system to reach zero roll 
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Position Feedback Controller. 

Position ftnd Velocity Feedback Controller. 
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angle position. Determination of Kl (proportional constant) 
was a trial and error process. Even though position feed- 
back alone gave big imorovement in the system dynamics, as is 
to be discussed in the following pages, it failed to stabil- 
ize the system in some operating conditions that are consider- 
ed very likely to be faced. This leads to the compensation 
of the system by velocity feedback and other modifications 
(Limiter in the position feedback channel). In the following 
pages the design procedure which uses extensive computer sim- 
ulation is discussed from the simplest case of proportional 
controller to the last modification that stabilized the system 
in various operating conditions and over a wide range of speeds. 

1. Proportional Controller 

The first attempt to control the roll was to design 
a proportional controller as is shown in Figure 57.1. Refer- 
ring to the figure, it is seen that ordered fairwater plane 
deflection is a function of the roll angle and the proportional 
constant Kl. Such that 

ORDERED FAIRWATER DEFLECTION = DFOD = -Kl*ROLL 
Since the reference signal is zero the system always looks for 
zero roll angle. Ordered fairwater deflection as an input to 
the fairwater planes actuator causes deflected planes at the 
output due to the actuator dynamics. Actuator dynamics were 
the same as those of the sternplanes actuator shown in Figure 
39. With the counter roll moment created differentially de- 
flected fairwater planes, system dynamics try to decrease the 
roll angle. 



96 



: 



Determination of Kl was a trial and error process. 

The stability range was 1 through 4. By inspection of the 
results of computer simulations, Kl = 3 was chosen as a best 
choice in the sense of snap roll and steady state value of 
the roll responses. With the determined Kl = 3 value the 
system was tested at 24 knots approach speed to a 35° rudder 
command and the results are shown in Figures 58 through 62. 
Before analyzing the results it is instructive to indicate 
here that all through the design procedure the worst condition 
that could happen was always taken into account. For this 
reason the controller was designed at 24 knots base speed and 
to a 35° constant rudder angle which could give the worst snap 
roll in the range of soeed of interest. By lookino at the 
Figures 58 through 62, the following results can be summarized: 
1. As was predicted, the snap roll decreased to almost 

4°. 6 from 37° and roll response reached a steady state 
value of 3°. 5. In Figures 52, 53, and 54 it was shown 
that before implementation of the roll controller the 
system was unstable to an identical test. Adding the roll 
controller made the system stable at high speed. In depth 
and pitch responses great improvement has been made. 

Steady state values of the depth and pitch were almost 
8 feet and 4°. 5 which were acceptable. While this depth, 
pitch and roll control improvement has been obtained, the 
stern and sailplane was used moderately such that they 
didn't reach any saturation (± 35°). 
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.Figure 58. 



Depth vs. Time. Kl = 3. UCK = 24 Knots. 

Rudder Ordered = 35°. Initial Roll Angle = 0°. 
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Figure 59. Pitch vs. Time. K1 = 3. UCK = 24 Knots. 

Rudder Ordered = 35°. Initial Roll Angle 
= 0 °. 
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Figure 60. Roll vs. Time. 

Ordered = 35°, 



K1 = 3. UCK = 24 Knots. Rudder 
Initial Roll Angle = 0°. 
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Figure 61. Sternplane Angle vs. Time. Kl = 3. UCK = 

24 Knots. Rudder Ordered = 35°. Initial 
Roll Angle = 0°. 
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Figure 62. Sailplane Angle vs. Time. Kl = 3. UCK = 
24 Knots. Rudder Ordered = 35°. Initial 
Roll Angle = 0°. 
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Figure 63, Depth vs. Time. Kl = 3. UCK = 24 Knots. 

Rudder Ordered = 35^. Initial ‘Roll Anale 
= -5°. 
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Figure 64. Pitch vs. Time. Kl = 3. UCK = 24 Knots. 

Rudder Ordered .= 35°. Initial Roll Angle 
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Figure 65. 



Roll vs. Time. Kl = 3. UCK = 24 Knots. 
Rudder Ordered = 35°. Initial Roll Angle 
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2. Investigation of the roll response in Figure 60 lead 

to further study. In the steady state value of the roll, 

the oscillation with almost fixed amplitude and fre- 

0 

quency was observed. The amplitude was around 0*1. As 
a first guess, the oscillation was thought to be a stable 
limitcycle. To expose the problem the system was tested 
with initial roll response (5^ inboard roll angle) at 24 
knots to a 35° constant rudder angle. Figures 63, 64, 
and 65 record the depth, pitch, and roll responses. The 
system became unstable. The system failed to compensate 
itself when it was commanded to turn with the 35° rudder 
deflection and 5° initial inboard roll. Any submarine 
maneuvering in the tactical area may experience a turn 
with initial roll angle. Disturbances from the heavy sea 
state might cause the submarine to roll eventhough it is 
in a straight course. If the submarine is commanded to 
turn at this moment, it was shown that its control is lost. 
That is why the situation was thought unacceptable and com- 
pensation of the system was attempted. It was thought that 
the reason for the failure was lack of enough feedback in- 
formation and compensation of the system with velocity 
feedback was attempted. 

The following section describes the modification and improve- 
ment obtained. 

2 . Compensation Of The System With Velocity Feedback 
The modified controller block diagram is shown in 
Figure 57.2., and DFOD = -K1*R0LL-K2*P . 
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It was felt that the velocity feedback should give to the 
system better dampina and better stability characteristics. 
Determination of the K2 value was again a trial and error 
process and the stability range was found to be 1 to 15 by 
computer simulation. By inspection of the results of the 
tests with different K2 values, K2 = 10 was considered the 
best choice in the sense of snap roll (max roll) and steady 
state value of the roll responses. Since the reason for intro- 
ducing velocity feedback was instability in the presence 
of an initial roll angle, tests with exaggerated initial roll 
(20° inboard) v/ere made. It was previously denoted that the 
ship has almost 37° snap roll at 24 knots to a 35° rudder 
command. To force the ship to turn with a hard rudder command 
when it already has a very serious initial roll angle was 
thought a good example of the capability to control the ship 
in three dimensions. In fact, any submarine cruising in a 
straight course is unlikely to face this much roll angle from 
heavy seastate. In Figure 66 the result of this test is shown. 
Before the implementation of the velocity feedback the system 
was unstable to a turn command with initial 5° roll angle. 

And now it is stable even with 20^ initial inboard roll. The 
result was such that steady state roll angle was almost 4°, 5 
and before reaching steady value on oscillation v;ith decreasing 
amplitude was observed. In all of the aforementioned tests, 
the ship response was investigated to a steady 35° rudder 
command. The rudder was commanded to a 35° deflection at the 
beginning of the simulation and after reachina 35° it was not 
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Figure 66. Roll vs. Time. K1 = 3, K2 = 10. UCK = 24 Knots. 

Rudder Ordered = 35°. Initial Poll Angle = -20°. 
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changed to any other command. But any submarine maneuvering 

in tactical areas can exoerience successive rudder commands 

in opposite directions. For this reason the test, in which 

o 

the rudder deflection was commanded to 35 in the time inter- 
val of 20 to 80 seconds was required. Figures 67 through 71 
record the depth, pitch, roll, sternolane and sailolane de- 
flection, In this test the initial roll angle was again 20° 
inboard. The result was unstable. Inspection of these 
figures reveals a very important reason of the failure. Curve 
number 2 in these figures represents the rudder response. But 
since in the vertical axis the automatic scaling was used 
associated with the output responses of interest, in some 
figures the rudder response in the time interval of 20 to 80 
seconds was seen less than 35°. For this reason, to overcome 
misinterpretation, curve number 2 should be interpreted as the 
time interval v;here rudder deflection was 35°, By inspection 
of the figures 67 and 68 its seen that the depth and oitch re- 
sponse are almost unchaged. The very small changes are due to 
the roll oscilation which stems from the initial roll angle. 

In Figure 70 it is shown that in the time interval of 0 - 20 
seconds the sternplanes oscillate with very small amplitudes 
to compensate these depth and pitch changes. But in the same 
time interval (0 - 20 seconds) the sailplanes oscillate with 
big amplitudes to overcome the roll oscillation started with 
the initial 20° inboard roll angle (Figure 71). Because of 
the roll controller, the initial roll angle has more effect 
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Figure 67.1, 

. 2 . 



Roll VS. Time. Kl = 3, K2 = 10. UCK = 24 
Knots. Initial Roll Angle = -20°. 

Rudder Response vs. Time. (Rudder Ordered = 35°). 
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Figure 68.1. Pitch vs. Time. Kl = 3, K2 = 10. UCK = 

24 Knots. Initial Roll Angle = -20°. 

.2. Rudder Response vs. Time. (Rudder Ordered 
= 35°). 
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Figure 69,1. Roll vs. Time. Kl = 3, K2 = 10. UCK = 
24 Knots, Initial Roll Angle = -20°. 

.2. Rudder Response vs. Time (Rudder Ordered 
= 35°. ) 
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40.00 80.00 120.00 160.00 200.00 

70,1, Sternplane Angle vs. Time, Kl = 3, K2 = 10, 
UCK = 24 Knots, Initial Roll Angle = -20°, 
,2, Rudder Response vs. Time (Rudder Ordered = 
35°) . 
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Figure 71.1. Sailplane Angle vs. Time. Kl = 3, K2 = 10. 

UCK = 24 Knots. Initial Roll Angle = -20°. 
.2. Rudder Resoonse vs. Time (Rudder Ordered = 
35°) . 
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on the sailplanes than on the sternplanes. In the tine inter- 
val of 0 - 20 seconds of Figure 69 the roll oscillation started 
with the initial roll angle is seen. Since the direct effect 
of the roll angle on the sailplane deflection, this initial 
roll oscillation causes the sailplane to oscillate and this 
oscillation causes initial roll oscillation not to be damned 
out. As this is going on, at time equal to 20 seconds the 
rudder was commanded to 35*^ full deflection and this gives more 
unstability. After analyzing this and previous results the 
following conclusion has been reached. Since the roll angle 
was directly fedback via Kl, the sudden and big roll response 
changes cause the sailplane to oscillate with big amplitudes 
and this leads to unstability. To overcome this difficulty a 
limiter was placed in the roll error feedback channel with the 
magnitude of ±5 (The same reason had lead to the placing 
depth and pitch error limiter in the depth and pitch controller 
design as was discussed in Section IV. B.). After this last 
modification the system was tested under various conditions 
which are to be discussed below. The complete controller block 
diagram with the depth and pitch controller is shown in Figure 
72. By inspection, the simulation results of the last version 
of the controller can be sumarized as follows: 

1. In Figures 73 through 87 the depth, pitch, roll, stern- 
plane and sailplanes responses at 24, 18, and 12 knots to 
35° rudder command in the time interval of 20 - 80 seconds) 
with the initial 20° inboard roll angle are shown. The 
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Figure 73.1. Depth vs. Time. Final Result With Roll 
Error Limiter. Kl = 3 , K2 = 10. UCK = 

24 Knots. Initial Roll Angle = -20°. 

.2. Rudder Resoonse vs. Time. (Rudder Ordered 
= 35°) . 
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Initial Roll Angle = -20°. 

.2. Rudder Response vs. Time (Rudder Ordered = 
35°) . 
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Figure 75.1. Roll vs. Time. Final Result With Roll Error 
Limiter. Kl = 3, K2 = 10. UCK = 24 Knots. 
Initial Roll Angle = -20°. 

.2. Rudder Response vs. Time. (Rudder Ordered 
= 35°) 
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Figure 76.1. Sternplane Angle vs. Time. Final Result With Roll 
Error Limiter. Kl = 3, K2 = 10. UCK = 24 Knots. 
Initial Roll Angle = -20°. 

.2. Rudder Response vs. Time (Rudder Ordered = 35°). 
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Figure 77.1. Sailplane Angle vs. Time. Final Result With Roll 
Error Limiter. Kl = 3, K2 = 10. UCK = 24 Knots. 
Initial Roll Angle = -20°. 

.2. Rudder Response vs. Time (Rudder Ordered = 35°). 
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Kl = 3, K2 = 10. UCK = 18 Knots. Initial '^oll 
Angle = -20°. 

.2. Rudder Response vs. Time (Rudder Ordered = 35°) 
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Figure 79.1. Pitch vs. Time. Final Result With Roll Error 
Limiter. K1 = 3, K2 = 10. UCK = 18 Knots. 
Initial Roll Angle = -20'^. 

.2. Rudder Response vs. Time (Rudder Ordered = 35°). 
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Figure 80,1. Poll vs. Time, Final Result With Roll Error 
Limiter. K1 = 3, K2 = 10, UCK = 18 Knots. 
Initial Roll Angle = -20°. 

,2. Rudder Response vs. Time (Rudder Ordered = 35°). 
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Figure 81.1. Sternplane Angle vs. Time. Final Result N.i th 
Roll Error Limiter. Kl = 3, K2 = 10. UCK = 

18 Knots. Initial Roll Angle = -20°. 

.2. Rudder Response vs. Time (Rudder Ordered = 35°). 
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Figure 82.1. Sailplane Angle vs. Time. Final Result With Roll 
Error Limiter. Kl = 3, K2 = 10. UCK = 18 Knots. 
Initial Roll Angle = -20°. 

.2. Rudder Response vs. Time (Rudder Ordered = 35°). 
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Figure 83.1. Depth vs. Time. Final Result With Roll Error 
Limiter. Kl = 3, K2 = 10. UCK = 12 Knots. 
Initial Roll Angle = -20°. 

.2. Rudder Response vs. Time (Rudder Ordered = 35°) 
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Figure 84.1. Pitch vs. Time. Final Fesult With Roll Error Limiter. 

K1 = 3, K2 = 10, UCK = 12 Knots. Initial Roll Ancle 
= - 20 °. 

.2. Rudder Response vs. Time (Rudder Ordered = 35°). 
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Figure 85.1. Roll vs. Time. Final Result With Roll Error Limiter 
K1 = 3, K2 = 10. UCK = 12 Knots. Initial Poll 
Angle = -20°. 

.2. Rudder Response vs. Time (Rudder Ordered = 35°). 



129 



i 

\ 



i 

\ 



35 ° 




Figure 86.1. Sternplane Angle vs. Time. Final Result With 
Roll Error Limiter. Kl = 3, K2 = 10. UCK = 

12 Knots. Initial Roll Angle = -20°. 

.2. Rudder Resnonse vs. Time (Rudder Ordered = 35°). 
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Figure 87.1. Sailplane vs. Time. Final Result With Poll Error 
Limiter. K1 = 3/ K2 = 10. UCK = 12 Knots. 
Initial Poll Angle = -20°. 

.2. Rudder Response vs. Time (Rudder Ordered = 35°). 
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Figure 88. Depth vs. Time. Final Pvesult With Roll Error 
Limiter. Kl = 3, K2 = 10. UCK = 24 Knots. 
PvUdder Ordered = 35°. Initial Roll Angle = 0°. 
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Figure 89. Pitch vs. Time. Final Result With Roll Error 
Limiter. K1 = 3, K2 = 10. UCK = 24 Knots. 
Rudder Ordered = 35*^. Initial Roll Anale = 0*^. 
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Figure 90. Roll vs. Time. Final Result With Roll Error 
Limiter. Kl = 3, K2 = 10. UCK = 24 Knots. 
Rudder Ordered = 35°. Initial Roll Angle = 0°. 
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Figure 91. Sternplane Angle vs. Time. Final 
Roll Error Limiter. Kl = 3, K2 = 
24 Knots. Rudder Ordered = 35°. 
Angle = 0°. 
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Figure 92. Sailplane Angle vs. Time. Final Fesult With 
Roll Error Limiter. K1 = 3, K2 = 10. UCK = 
24 Knots. Rudder Ordered = 35*^. Initial 
Roll Angle = 0°. 
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Figure 93. Depth vs. Time. Final Result With Roll Error 
Limiter. Kl = 3, K2 = 10. UCK = 18 Knots. 
Rudder Ordered = 35°. Initial Roll Ancle = 0°. 
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Figire 94. Pitch vs. Time. Final Result Nith Roll 
Error Limiter. Kl = 3, K2 = 10. UCK = 

18 Knots. Rudder Ordered = 35^. Initial 
Roll Angle = 0°. 
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Figure 95. Roll vs. Time. Final Result With Roll Error 
Limiter. Kl = 3, K2 = 10. UCK = 18 Knots. 
Rudder Ordered = 35^. Initial Roll Angle = 0*^. 
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Figure 96. Sternplane Angle vs. Time. Final 
Roll Error Limiter. Kl = 3, K2 = 
18 Knots. Rudder Ordered = 35°. 
Angle = 0°. 
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Figure 97. Sailplane Angle vs. Time. Final Result With 
Roll Error Limiter. Kl = 3, K2 = 10. UCK =. 
18 Knots. Rudder Ordered = 35°. Initial 
Roll Angle = 0°. 
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Figure 98. Depth vs. Time. Final Pesult With 

Roll Error Limiter. K1 = 3, K2 = 10. UCK = 

12 Knots. Rudder Ordered = 35°. Initial Roll 
Angle = 0°. 
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Figure 99 . Pitch vs. Time. Final Result With Poll Error 
Limiter. Kl = 3, K2 = 10. UCK = 12 Knots. 
Rudder Ordered = 35°. Initial Roll Angle = 0°. 
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Figure 100. Roll vs. Time. Final Result With Roll Error 
Limiter. K1 = 3, K2 = 10. UCK = 12 Knots. 
Rudder Ordered = 35°. Initial Roll Angle = 0°. 
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Figure 101. Sternplane Angle vs. Time. Final Result With 
Roll Error Limiter. Kl = 3, K2 = 10. UCK = 

12 Knots. Rudder Ordered = 35*^. Initial Poll 
Angle = 0°. 
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Figure 102. Sailplane Anqle vs. Time. Final Result With 
Roll Error Limiter. Kl = 3, K2 = 10. UCK = 

12 Knots. Rudder Ordered = 35°. Initial i^oll 
Angle = 0° . 
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Figure 103. Roll vs. Tiine. Final Result With Roll Error 
Limiter. Kl = 3, K2 = 10. UCK = 24 Knots. 
Rudder Ordered = 35°. Initial Roll Anale = 
- 20 °. 
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Figure 104. Roll vs. Time. Final Result With Roll Error 
Limiter. Kl = 3, K2 = 10. UCK = 18 Knots. 
Rudder Ordered = 35°. Initial Roll Angle 
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105. Roll vs. Tiitie. Final Result With Poll Error 
Limiter. K1 = 3, K2 = 10. UCK = 12 Knots. 
Rudder Ordered = 35°. Initial Roll Ancle = 
- 20 °. 
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same test without roll limiter as was discussed above 
gave an unstable system but after inserting the limiter 
the result was stable over the speed range of interest. 

In contrast to the results shown in Figures 67 through 71, 
there is no oscillation in the roll and sailplanes re- 
sponses in the time interval of 0 - 20 seconds and it 
was observed that the limiter gave very good damning to 
the system. This last statement can be justified by 
inspection of Figure 75 such that at time equal 20 
(when the rudder is commanded to full deflection) the 
roll response is almost zero without any previous oscil- 
lation. As the stabilization in three dimensions v/as being 
reached the stern and sailplanes never went into saturation. 

2. Figures 88 through 102 record the depth, pitch, roll, 
sternplane and sailplanes responses at 24, 18, and 12 
knots to 35° rudder command (the rudder is commanded to 
full deflection of the beginning of the simulation) with 
zero initial roll angle. It is observed that without 
using excessive stern and sailplane deflection great sta- 
bility in the three dimension has been reached. As was 
shown before. Figures 44 through 55 record the responses 
of the system without roll controller to the identical 
test. Comparison of these tests is made in the Table II. 

It is seen that at high speed great improvement in the 
sense of maximum and steady state value has been obtained. 
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TABLE II 



IMPROVEMENT DUE TO ROLL CONTROLLER 



UCK 




WITHOUT ROLL CONTROLLER 


WITH ROLL CONTROLLER 


MAX 


STEADY STATE 


MJIlX 


STEADY STATE 




Depth Changed 






8.5 ft. 


8 ft. 


24 


Pitch 




UNSTABLE 


5°. 5 


4°. 4 




Roll 






5°. 3 


4° 




Depth Changed 


10.1 


8.7 


8.5 


8 ft. 


18 


Pitch 


9.5 


5° 


5°. 6 


4°. 7 




Roll 


23° 


6° 


4°. 6 


3°. 2 




Depth Changed 


9.7 ft. 


9.5 ft. 


9.2 


9 


12 


Pitch 


7° 


5°. 7 


6.7 


5.2 




Roll 


10°. 3 


2.5 


3°. 5 


2° 
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3. Figures 103 through 105 record the roll resoonses of the 
system at 24, 18, and 12 knots to 35° rudder deflection 
(the turn was commanded at the beginning of the simu- 
lation) v'ith 20° inboard roll angle. The results are 
stable. Comparison of Figure 103 with Figure 66 shows 
the improvement obtained by inserting a roll limiter such 
that it damped out the oscillation. 

It was shoT^m that utilizing the fairwater planes as 
the control surface of the roll controller not only stabilized 
the roll response but also gave big improvement in the depth 
and pitch controller. In the following section stability 
tests of the system to disturbance moments is discussed. 

C. ROLL CONTROLLER STABILITY TESTS 

Since the roll controller design was carried out usinc 
nonlinear equations in six degrees of freedom, the stability 
analysis of the system v;as not practical. To investigate the 
stability comouter simulation was used by applying disturbance 
moments. The magnitude of the moments was 
DISTURBAxNCE = KS*(UCK)^*L^ where 
KS = 4.0*0.0003*576/(UCK)^ 

The disturbance moment was calculated by assigning non- 
zero value of KS which was normally zero and reoresents the 
hydrodynamic coefficient which gives the rolling moment when 
body angle and control surface angles are zero. The magnitude 
of the disturbance was four times bigger than that used in the 
stability analysis of Stamps controller. At 24 knots the 
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magnitude of the disturbance was 1.1028x10^ ft-lb. such that 
the Stamps controller was unable to stabilize the system to 
the disturbances of the magnitude in the order of 10^. Figures 
106 through 120 record the depth, oitch, roll, sternplane, and 
sailplanes responses to the inboard roll moment dis- 
turbance mentioned above. The distrubance was applied at time 
equal 20 as a step moment. Curve number 2 represents the dis- 
turbances. Inspection of the results shows that the system is 
stable in three dimensions to the disturbances moment. The 
worst affect of the disturbances was seen at low speed test 
v/hich was normal. But even in this case, after a transient re- 
sponse the system was stable. 
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Figure 106.1. 

. 2 . 



Depth vs. Time. Roll Stability Steo Test. 

UCK = 24 Knots. Rudder Ordered = 35*^. 
Disturbance Moment vs. Time (1. I028*10^ft-lb . ) . 
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Figure 107.1. 

. 2 . 



Pitch vs. Time. Roll Stability Steo Test. 
UCK = 24 Knots. Rudder Ordered = 35°. 
Disturbance vs. Time. (1. 1028*10^ft-lb. ) . 
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Figure 108.1. 

. 2 . 



Roll vs. Time. Roll Stability Step Test. UCK 
= 24 Knots. Rudder Ordered = 35°. 

Disturbance Moment vs. Time (1. 1028xl0^ft-lb. ) . • 
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Figure 109.1. 

. 2 . 



Sternplane Ancle vs. Time. Roll Stability Steo 
Test. UCK = 24 Knots. Rudder Ordered = 35°. 
Disturbance Moment vs. Time (1. 1028xl0^ft-lb. ) . 
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Figure 110.1, 



. 2 . 



Sailplane Angle vs. Time. 
Test. UCK = 24 Knots. 
Disturbance Moments vs. Time. 



Poll Stability Step 
Rudder Ordered = 35°. 



(1. 1028xlo'ft-lb. ) 
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Figure 111.1. Depth vs. Time. Roll Stability Step Test. 

UCK = 18 Knots. Rudder Ordered = 35°. 

.2. Disturbance Moment vs. Time. (1. 1028xl0^ft-lb. )■. 
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Figure 112.1. 
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Pitch vs. Time. Roll Stability Steo Test. 

UCK = 18 Knots. Rudder Ordered = 35°. 

Disturbance Moment vs. Time. (1. 1023xlo"^ft-lb. ) . 



160 



10.00 - 8.00 - 6.00 - 4.00 - 2.00 0,00 2.00 4.00 



X SCALE — 40 Seconds Per Inch. 

Y SCALE (1) =2 Degrees Per Inch 




■/ 



1 



0.00 



TOTW 



80 ": 00 



120 . 00 ^ 160.00 



200.00 



Figure 113.1. Roll vs. Time. Roll Stability Step Test. 

UCK = 18 Knots. Rudder Ordered = 35^. 

.2. Disturbance Moment vs. Time. (1. 1028xlo'^ft-lb. ). 
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Pigure 114.1. 

. 2 . 



Sternplane Angle vs. Time. Roll Stability Steo 
Test. UCK = 18 Knots. Rudder Ordered = 35°. 
Disturbance Moment vs. Time. ( 1. 1028xlo'^f t-lb . ). . 
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Fiaure 115.1. Sailplane Angle vs. Time. Roll Stability Steo 

Test. UCK = 24 Knots. Rudder Ordered = 35°. 

.2. Disturbance Moment vs. Time. (1. 1028xlo"^ft-lb. ) . 
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Figure 116.1. 



Depth vs. Time. 
UCK = 12 Knots. 



Roll Stability, Step Test. 
Rudder Ordered = 35°. 



.2. Disturbance Moment vs. Time. (1. 1028xlo'^ft-lb. ) 
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Figure 117.1. 
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Pitch vs. Time. Roll Stability Steo Test. 

UCK = 12 Knots. Rudder Ordered = 35*^. 

Disturbance Moment vs. Time. (1. 1028xlo'^ft“lb. ) . 
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Roll vs. Time. Roll Stability Sten Test. 

UCK = 12 Knots. Rudder Ordered = 35°. 

7 

Disturbance Moment vs. Time. (1.1028x10 ft-lb.). 
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Figure 119.1. 
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Sternplane Angle vs. Time. Roll Stabilitv 
Step Test. UCK = 6 Knots. Rudder Ordered = 35^ 

7 

Disturbance Moment vs. Time. (1.1028x10 ft-lb.) 
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Figure 120.1. 
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Sailplane Angle vs. Time. Roll Stability Sten 
Test. UCK = 12 Knots. Rudder Ordered = 35°. 

7 

Disturbance Moment vs. Time. (1.1028x10 ft-lb.) 
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VI. DISCUSSION, CONCLUSIONS AND RECOMMENDATIONS 

FOR FURTHER WORK 



A. DISCUSSION AND CONCLUSION 

The characteristic of a submarine which is totally 
submerged and in a high speed turn was discussed in Chapter II 
and the simulation results were given in Chapter III. The 
unhydrodynamic body structure when it is submerged, v;hich is 
mainly due to the appendages (sail) , was considered the main 
source of the problem. The necessity of the sail that provides 
rooms for periscope and other vital controls that are located 
in the tower does not oermit a small sail structure. On this 
account, a control scheme must be used to compensate the prob- 
lems that stem from the sail. The complexity of controlling a 
submarine in a high speed turn, comes from the coupling between 
states and three-dimension dynamics. In Chapter IV an automatic 
depth and pitch controller was presented briefly with the modi- 
fication of the original design represented in Reference 7. 

The automatic control scheme gave the advantage of taking depth 
and pitch error and their rate into account. But, since only 
the sternplane was used as a control surface, the results made 
the roll controller design necessary for controlling the ship 
within the wide range of speeds and under various disturbance 
conditions. The imposition of the fairwater plane to control 
roll, which was the main goal of this thesis, gave a great 
amount of expected improvement to the system. The roll 
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controller scheme was carried out by using equations of motion 
(nonlinear) in ^<;hich the vertical plane dynamics were stabil- 
ized by the imposition of the depth-pitch controller, and roll 
state and roll rate gave enough information to provide closed- 
loop control for the ship. The control scheme was a direct 
roll controller compared to the controller designed by Stamos 
in Reference 2. Stamps's controller was an indirect roll con- 
troller and in his design the rudder was used as a function of 
the integral of roll error which was the difference between the 
allowed maximum roll and the actual instantaneous roll angles. 
The initial rudder ordered was to be such that the peak roll 
expected for a given approach speed would be less than the 
maximum allowed roll. The integral of roll error v;as then to 
be computed and scaled to represent an additive tern applied 
to the initial rudder order (9ro). The instantaneous rudder 
order (DRCOM) was then to be computed 

DRCOM = 6ro + K_ ((^ -d) dt. 

I r max r act 

The technique reported in this thesis has the benefits of 
simplicity and costs much less than alterations of the design 
criteria of Navy submarines. But since it prohibited use of 
a hard rudder deflection (35^) in the early phase of turn, the 
time required to pass any specified yaw angle v;as quite long. 

For this reason, the roll control, in a direct sense, by making 
use of fairwater planes as control fins, w'as investigated in 
the proposed design. As a result, not only were the dangerous 
snap rolls decreased to magnitudes of 5°, but also with the hard 
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rudder deflection a big yaw rate was achieved. In Table III, 
the performance of the submarine with each of these two roll 
controllers is compared. In the Stamps* design, maximum roll 
angle is predetermined and the controller is designed not to 
exceed this value. The higher allowed maximum roll angle gives 
higher yaw rates. In the proposed desian, maximum roll angle 
is the result of the system's dynamics. 

TABLE III 

COMPARISON OF THE PROPOSED AND STAiMPS ' DESIGN 



TYPE OF 


MAX 


APPROACH 


TIME REQUIRED IN 


1 

SECONDS TO 


CONTROLLER 


ROLL 


SPEED 


PASS SPECIFIED YAW ANGLE 








45° 


90° 


180° 


STAMPS* DESIGN 




24 


52 


91 


7100 




7°. 5 


18 


39 


67 


7100 






12 


27 


47 


83 


DIFFERENTIAL 


5°. 61 


24 


11.5 


20 


41 


FAIRWATER 


4°. 69 


18 


14 


25 


54 


PLA1<IE 


3°. 4 


12 


19 


36 


80 



It is seen that the use of differential fairwater planes 
provides a faster response with reduced maximum roll angle at 
all speeds. 

In the proposed design, only one fairwater plane actu- 
ator was used and it was assumed that it gives deflection in a 
differential mode. Switching criteria from the conventional 
usage of the fair^^ater plane (both starboard and port sides 
move simultaneously in the same direction) to the differential 
mode was not discussed. 
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This research has shov/n that imposition of the faiIn^7ater 
olane as a part of the roll controller improved the control of 
a turning high speed submarine in three dimensions. Dangerous 
snap roll was decreased to around 5° and with only sternplane 
and the stabilized roll response high improvement in the depth 
and pitch characteristic was achieved. 

B. RECOMMENDATIONS FOR FURTHER WORK 

The follo^ving research xvould be worthwhile in future 
studies : 

1. Automatic deoth and pitch controller by using only 
sternplane in the sense of an optimal feedback controller 
for linear tracking should be designed and the combina- 
tion of this design with the present roll controller 
should be investigated. This study would give good in- 
sight into the problem which comes from the depth changing 
when the submarine is in a high speed turn. 

2. Switching criteria from the roll controller of Reference 
2 to the proposed roll controller should be established 
and these two designs should combine together. This would 
give better depth and oitch control to the ship when it 
has low speed where a great amount of roll control is 
not needed. 

3. With the coordination of NSRDC a more accurate hydrodynamic 
coefficient associated with the differentially deflected 
sailplane should be obtained. 
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4. Addition of an integral of error to the roll control 
should be studied to see if the steady state roll angle 
can be reduced. 
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. APPENDIX A 



EQUATIONS OF MOTION 



The followlns 5 set of equations are referred to a body 
fixed systen of axes which arc coincident with the prin- 
cipal axes of Inertia of the body. The origin of this 
axis-system Is located at the assumed center of mass of 
the bodv 



Equation of Motion Along Che Body Axis System x-Axis 



n,(u - vr ► wq). I I* [x ♦ X,^ ' ■■■ ' ’t] 




t |iV [x 




+ 



vvn* 



’ (n' - 1) V* 




' (n» - I ) w* 



*|<vx 



6s6sn' 



• in' - 1) 




' (n* - 1) 



sin 0 
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Equation of Motion Along the Body Axis System y-Axls 



m(v - wp > ur) - 

* I'* [^p, ■P’’* ’^pIpI 'pIpI] 

[r. '..ptY^i^i.iXiKv'tw’iil |,|] 

♦ |J» [y, '“p^^|p|6p '-IrlirtYj, ’up] 

'^Tu- ' 

+ -|j‘ fY , ' u‘ <• Y ^ ' uv t Y ^ 1^1 ' v|(v‘ f w'l^l] 

' ‘r 

♦I''"' ^ 6r«' ‘r 

♦I'' '■.n' '<■*■- "'“Y 

*Z‘‘ '^wv’y'y' 

+ -f X* (FJ V* + w* (-w) Sin 0)t f Multiplied by 

^ y vs 



^ EW j sin <p cos 0 



for large angles of 
attack near -90® 
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Equaclon of Motion Along the Body Axis System z-Axl# 



m(w - uq + vp) = 2 ^ ^ 4 



I'* [^„ ’'P'] 



Note L 



w ’ ^ ^ vr ' ^ ^ vp ' vp ' 



+ 1 ^ 



* [z* '“*+Zw ’“'^^^wlwl +'*'*)^l] 

+ 1^* f^|wj’ “!'"! " ^ww’ w^)^[ + V*#] 

+ |ZV [Z^, Ms.Z^^ -6b] 

^ 1^* [^ wa' ' f"' - l)uw+ Z • (n« - 1) wl(v» + w*)^!] 






+ -?i* (F ) V* + w’ 

2 ' x'vs JJ— 

+ EW. cos 9 cos 4> 



V Sin cot 



# Multiplied by 

u 

u 

for large angles of 
attack near -90® 

Note 1 

when not multiplied by 
u add to Z ' 

T3 vp 
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Equation of Motion About the Bod/ Axis S/scern x-Axls 



IjjP + (I* - ly) = I ^ ' P ^ ^ qr ' ^ * 



+ I* Fk 'up+K 'ur + K. 'v + 

6 I* P * ^ 



+ -| [k , ' u* + K y • uv + K ^1 • V 



+ •? K ' vw 

Z VW 



♦ 1^’ “‘K»r ■ K 



+ Bz Q ein <{) coa 8 



^ pIpI ' pIpI] 

l(v* 
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Equation of Motion About the Body Axta System y-Axls 



Note 1 

lyq + (I, - rp = [m . • q + M • r^l f M ‘ rp + AM • rp#] 



+ -| A* [m ^ ' uq + M 1 ^ 1 ^^ • u|ql6s + M ' 1< V* + 

•f Z* ^ *w+M^^ 'vrrf+M^p 

♦ l^*^qn' ’ 

+ |X» [m, *u*+M^ 'uw+M^,^, • wl{v»+w*)^r] 

+ -|x^ [M|^| ’ ul w( + |w (v* + w^)^| f V*? 1 

"I"' [^6« ’‘a^^6b *«b] 



+ 



-fx’ M , 
2 wn* 



(n* • 1) uw 



* l'“ “w|wln’ ’ '"' - ‘> 



+ Bx sin 0 



EW ^ x^. cos 9 cos ^ 



# Multiply by g for 
large ancles ofaltack 
near -90^ 

Note 1 

when not multiplied 
by ii add to M * 

U 
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» 1 * 



HquaClon of Hoclon AbouC Che Bod/ Axis S/scea t-Axls 



+ (I 



y 



•'x'PI [«} 'p] 

+ !<• [n^ '"P*”|r| 6 , ' N|,|, ' “'ihr] 

* 2 '“P*Mj ' '^1 

*T'* ”rn' ' 

+ [n. 'u^+N^ 'uv+N^i^i • v1(v» + I] 

+ • (ti* . 1 ) uv 

* ^v|vln' ' vl(v* +w*)h 

+ N • wvf 

^ wv 

+ EW^ x^. cos 0 Sincj) f Multiply by ^ for large 

angles of attack near - 90 ^ 
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AUXILIARY EQUATIONS 



^ s p ^ i sin9 

6 * (q - ^COs9 sin4> ) / C03 4> 

i> s ( r + 0 sin 4> ) / cos 8co3<() 

Xo - u CO 80 CO sii) + v(s in <J)sin0 CO CO s 4^3 in 

+ w ( 8in<{) sin^ f cos<{) sin 0 cos^) 

Yq • ucos9sini() + v{cosc()cos^J) + 3in4)sin03inf{)) 
4 w(cos<})sin0 sin^ - sin(}) cos ^ ) 

Zq = - u sin 8 + V cos 9 sin<{) + w cos 0 cos ({> 

U s (u* + + w*)^ 



(Fjj)p s ^ u* ai' + ag’ n’ + a,' n'* J 



when < n* 



I [ V + ha* n* + ba* n’^ ] when < n’ < 

^ X^ ^ c^’ + C3* n’ + Ca* J when k^ < n* < 3^ 

^ X^ ^ di + dj* n* + dj* J when n* < ka 



• • • 

*»!• '»3 

• t , 

« 1 » « 2 » «3 

Iff 

di, d2f dj 



Sets of non-dlrenslor.nl cocrriclcnts »iscd In the pro- 
pulsion equation above. set vhlcb. will be In effect 

at any tl.TC durlnrj a sln ilaVci maneuver will depend on the 
value of n and U'.e nuT^bers k^,k2»^2* 
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nomemclaturr: 



All symbols used in Che equations of motion and In the 
auxiliary equations and relationships which appear in this 
report are defined below. Any dimensions involved will be 
consistent with the f oo t -pOMnd-sccond system of units. All 
angles are in degrees. The Fortran variables corresponding 
to these symbols arc shown in Appendix 3 



SYMBOL 



B 



Si 

u 



U 



DEFINITION 

A dot over any symbol signi- 
fies differentiation with res- 
pect to time. 

Buoyancy force which is posi- 
tive upwards. 

Mass of the submarine including 
the water In the free ^loo^^ng 
spaces « 

Overall length of the submarine 

Linear velocity of origin of 
body axes relative to an earth- 
fixed axis system. 

Component of U along the body 
x-axls • 



V Component of U along the body 

y-axls. 

« 

V Component of U along the body 

£-axls . 
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Commaad speed: A sceady value 

of u for a given propeller rpta 
whence, |S) and control surface 
angles are zero. Sign changes 
with propeller reversal. 

Longitudinal axis of the body 
fixed coordinate axis system. 

Transverse axis of the body 
fixed coordinate axis system. 

Vertical axis of the body fixed 
coordinate axis system. 

Distance along the x^ axis of ao 
earth-fixed axis system. 

Distance along the y^ axis of an 
earth-fixed axis system. 

Olgcance along the axis of an 
earth-fixed axis system. 

Component of angular velocity 
about the body fixed x-axls. 

Component of angular velocity 
about the body fixed y-axls* 

Component of angular velocity 
about the body fixed z-axls. 

The z^coordlnate of the center 
of buoyance (C3) of the subma- 
rine. 



o< 


Angle of attack. 
Angle of drift. 


Sb f D& 


Deflection of bouplane ( or 
•ellplane } 


Sr j Dr- 


Deflection of rudder. 


5s , 05 


Deflection of aternplane. 


• 

a 


The ratio u /u. 

c 


0 

$ 

<!> 

C 


Pitch angle. 

Taw angle. 

Roll angle. 

Mass density of sea water. 


“i 


Weight of water blown from a 
particular ballast tank Ident- 
ified hy the Integer assigned 
to the Index 1* 


Oi 


Angular velocity. 


t 


Time. 


*el 


hocatlon along the body x-axls 
of the center of mass of th# l ^h 

ballast tank when this tank la 
filled with sea water. 
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I 







I 



X 



Propulsion force (see auxl- 
Ilnry equations and relation* 
ships). 

Moment of Inertia of a sub* 
marine about the x*axts. 



X 



y 



Moment of Inertia of a sub- 
marine about the y-axls. 



I 



z 



Moment of Inertia of a sub- 
marine about the r-axls. 



^p|p| '^qr* 







/ 

v» 




* ^ Y- / 

^vv» 



^Sr 



Non-dimensional constants each 
of which is assigned to a parti- 
cular force term in the equation 
of motion about the body x-axis. 



M. • 



M 



rr 



M • 

rp 



AM • 

r? 



M 



I X t 

,q| As* 



"»'■ «vr’- "vp'' ",/• 



V W W W WW VV OS 



”6b’* ”wn*’* 



V I V n 



•Asn^ 



Non-d 1 mens 1 ona 1 constants each 
of which Is assigned to a parti- 
cular force term in the equation 
of motion about the body y-axis. 
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V 



N 



c’* '"r’’ ''|r|4r’* '*M r’ 



%’• ”vp*’ '‘rn-** 

^vIvT' ^4r’’ ^4rn*’* ^'vn’’’ ”v|v!n'’ 



Non*d Inena Iona I conscancs each 
of which la as/vlgncd to a parti* 
cular force tern In the equation 
of motion about the body z*axls* 



5 • 

vv 






t 



T • X • X * X ’X • Non-dimensional constants each 
uu • vv • ww • 6 r4 r • 5 s6 a ’ 

of which la assigned to a parti- 

X-.,.’’, X •, X I ’ . X. . cular force term In the equation 

0 bo b vvn wwn da6an 

of notion along the body x-axls* 

*4r«rn-’ 



Y 



i 



» 

» 








t 



Mrl’’ %'• 






\n-'- ’'.Ma-'- <V 



VS 



Non-d Incns 1 ona 1 constants each 
of which is assigned to a parti* 
cular force term in the equatlo* 
of motion along the body y-axis 



2. • , Z * 
q rr 



^*rp’* ^w’* ^vr* ^vp'* 



A 

Z 

z 

z 



~vp * ^q * ^1 q| 4 s’ ^w| ql ’ 





. z.’ , z ’ , 


Z 1 1 ’ . 




qn* 


• a ' V 


w| wj 


I«I 


f 

vv 


2 • Z • 

’ *vv ’ 4s 


• ^4b’’ 


Z t 
vn 


v| V 


Jn” ^4sn” 













Non-dimensional constants each 
of which is assigned to a parti- 
cular forc«« term in the equatlo* 
of motion along the body z-axls 
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